Electron Transport in Silicon Nanowires: 
The Role of Acoustic Phonon Confinement and Surface Roughness Scattering 
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We investigate the effects of electron and acoustic-phonon confinement on the low-field electron 
mobility of thin square silicon nanowires (SiNWs) that are surrounded by Si02 and gated. We 
employ a self-consistent Poisson-Schrodinger-Monte Carlo solver that accounts for scattering 
due to acoustic phonons (confined and bulk), intervalley phonons, and the Si/Si02 surface 
roughness. The wires considered have cross sections between 3x3 nm 2 and 8x8 nm 2 . For 
larger wires, as expected, the dependence of the mobility on the transverse field from the gate 
is pronounced. At low transverse fields, where phonon scattering dominates, scattering from 
confined acoustic phonons results in about a 10% decrease of the mobility with respect to the bulk 
phonon approximation. As the wire cross-section decreases, the electron mobility drops because the 
detrimental increase in both electron-acoustic phonon and electron-surface roughness scattering 
rates overshadows the beneficial volume inversion and subband modulation. For wires thinner 
than 5x5 nm 2 , surface roughness scattering dominates regardless of the transverse field ap- 
plied and leads to a monotonic decrease of the electron mobility with decreasing SiNWs cross section. 
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I. INTRODUCTION 



Among the emerging devices for the future tech- 
nology nodes, silicon nanowires (SiNWs) have attracted 
much attention among researchers due to their potential 
to function as thermoelectric coolers?^ logic devices^ 
interconnects^ photodetectorSf 7 - as well as biological 
and chemical sensors. Although considerable work has 
recently been done on SiNW^i 9 -^ there is no consensus 
on these structures' electronic properties. 

The low- field electron mobility is one of the most im- 
portant parameters that determines the performance of 
field-effect transistors (FETs), thermoelectric (TE) cool- 
ers and sensors. Although its importance in ultra-short 
channel MOSFETs has been debated ) 11 ' 12 it certainly af- 
fects the conductivity of the wire interconnects, the figure 
of merit of TE coolers and the responsiveness of nanowire 
sensors. The study of the electron mobility in SiNWs so 
far has been inconclusive: Kotlyar et alM- and Jin et aW* 
have shown that the mobility in a cylindrical SiNW de- 
creases with decreasing diameter, whereas the works of 
Sakaki, 13 Cui et al-r- and Koo et alM- show higher mobil- 
ity in SiNWs compared to bulk MOSFETs. The contra- 
diction stems from two opposing effects that determine 
the electron mobility as we move from 2D to ID struc- 
tures: one is a decrease in the density of states (DOS) 
for scattering^ 3 - that results in reduced scattering rates 
and thereby an enhancement in the mobility; the second 
is an increase in the so-called electron-phonon wavefunc- 
tion overlap^ that results in increased electron-phonon 
scattering rates and consequently lower mobility. While 
important, these two competing phenomena do not paint 
a full picture of low- field transport in SiNWs, in which 
the effect of spatial confinement on the scattering due 



to surface roughness and acoustic phonons must be ad- 
dressed. 

Surface roughness scattering (SRS) is by far the 
most important cause of mobility degradation in conven- 
tional MOSFETs at high transverse fields. Intuitively, 
one would expect the SRS to be even more detrimental 
in SiNWs than conventional MOSFETs because SiNWs 
have four Si-Si02 interfaces, as opposed to one such in- 
terface in conventional MOSFETs. Although recent work 
has shown that the SRS is in fact less important in SiNWs 
than in bulk MOSFETs due to a reduction in the DOS^ 
and the onset of volume inversion (redistribution of elec- 
trons throughout the silicon channel)?^ a detailed study 
of the SRS-limited mobility in cylindrical SiNWs by Jin 
et alr^ shows a rapid monotonic decrease of mobility for 
wire diameters smaller than 5 nm. 

Scattering due to acoustic phonons is significantly 
altered in nanostructures due to the modification of the 
acoustic phonon spectrum in them. Extensive work 
on the effects of acoustic phonon confinement in III-V 
based nanostructure a 17 ' 18 shows a lower acoustic phonon 
group velocity ) 19 ' 20 lower thermal conductivity ) 21 ' 22 ' 23 
and increased acoustic phonon scattering rate a 24 ' 25 in 
nanoscale devices compared to their bulk counter- 
parts. As for silicon nanostructures, there is ex- 
perimental evidence of the acoustic phonons confine- 
ment in nanomembranesj2£ while recent works on SOI 
MOSFETs^ 7 -^^ and SiNWs^ . also show that the 
bulk phonon (linear dispersion) approximation underes- 
timates the scattering rates. 

In this work, we calculate the electron mobil- 
ity of gated square SiNWs by using a self-consistent 
Schrodinger-Poisson-Monte Carlo simulator and account- 
ing for electron scattering due to acoustic phonons (con- 



fined and bulk), intervalley phonons, and imperfections 
at the Si/SiC>2 interface. The wires considered have cross 
sections between 3x3 nm 2 and 8x8 nm 2 . Bulk-silicon 
effective mass parameters are used in the calculation of 
the scattering rates. The confined acoustic phonon spec- 
trum is obtained by using the xyz algorith m 31 ! 32 and 
the unscreened SRS is modeled using modified Ando's 
model^ 3 . that accounts for the finite thickness of the sili- 
con layer. For larger wires, as expected, the dependence 
of the mobility on the transverse field from the gate 
is pronounced; at low transverse fields, where phonon 
scattering dominates, scattering from confined acoustic 
phonons results in about a 10%-decrease of the mobil- 
ity with respect to the bulk phonon approximation. As 
the wire cross section decreases, the electron mobility 
drops because the detrimental increase in both electron- 
acoustic phonon and electron-surface roughness scatter- 
ing rates overshadows the beneficial effects of volume in- 
version (redistribution of electrons in real space) 16 i 34 and 
subband modulation (redistribution of electrons among 
different subbands ) 35 i 36 i 37 . For wires thinner than 5x5 
nm 2 , surface roughness scattering dominates regardless 
of the transverse field applied and leads to a monotonic 
decrease of the electron mobility with decreasing SiNWs 
cross-section. 

The electronic bandstructure in SiNWs is altered 
from that of bulk silicon due to the Brillouin zone 
folding i 30 i 38 i 39 The degeneracy of the conduction band 
minima in SiNWs is split; A4 valleys (four degenerate 
A valleys whose long axis is perpendicular to the SiNW 
axis) are found to be at the T point and the A2 val- 
leys (two degenerate A valleys whose long axis is par- 
allel to the SiNW) are at k x = ±0.37?r/a. The band- 
structure modification will certainly affect the interval- 
ley scattering rates in SiNWs, but deformation poten- 
tials and phonon energies needed to describe the inter- 
valley scattering are still only available for bulk silicon. 
Therefore, in order to be consistent, we do not employ 
the exact nanowire bandstructure, but rather rely on the 
bulk Si bandstucture and then solve the two-dimensional 
Schrodinger equation within the the envelope function 
and effective mass framework, with the bulk Si effective 
mass parameters. This approximation was proven ade- 
quate down to the 3 nm wire diameter i 3 ^ The splitting of 
the sixfold degenerate A valleys of bulk silicon into A2 
and A4 upon confinement is automatically accounted for 
by including the anisotropy of the electron effective mass 
(see Sec. UTTAl . 

This paper is organized as follows: Section [IT] de- 
scribes the device structure used in this study and the 
components of the simulator developed to calculate the 
mobility. A description of the acoustic phonon spectrum 
calculation in SiNWs is given towards the end of Sec- 
tion [TTJ In Section IIII1 we emphasize the importance 
of accounting for the acoustic phonon confinement when 
calculating the mobility in SiNWs and then present the 
results for the variation of mobility in square SiNWs with 
increasing spatial confinement. We conclude this paper 
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FIG. 1. (Color online) Schematic of the simulated 8x8 nm 2 
SiNW on ultrathin SOI. 



with a brief summary of the findings of our work (Sec- 
tion [IV| and a detailed derivation of the scattering rates 
due to acoustic phonons, intervalley phonons, and surface 
roughness (Appendices \M and [Bj . 



II. MOBILITY CALCULATION 
A. Device Structure and Simulator Components 

A schematic of the device considered in our study is 
shown in the top panel of Fig. [TJ It is a modified version 
of the ultra-thin, ultra-narrow SOI MOSFET that was 
originally proposed by Majima el a/4£ The thickness of 
the gate oxide, buried oxide, and bottom silicon substrate 
are 25 nm, 80 nm, and 700 nm, respectively. The trans- 
verse dimensions of the silicon channel are varied from 
8x8 nm 2 to 3 x3 nm 2 . For all the device cross sec- 
tions considered, the width of the oxide on both sides of 
the channel is 200 nm, the channel is doped to 3 x 10 15 
cm -3 , and the channel is assumed to be homogeneous 
and infinitely long. 

The simulator developed to calculate the electron 
mobility has two components: the first is a self-consistent 
2D Poisson-2D Schrodinger solver and the second is a 
Monte Carlo transport kernel. The former is used to cal- 
culate the electronic states and the self-consistent poten- 
tial distribution along the cross section of the wire and 
the latter simulates the transport along the wire axis. 
The finite barrier at the Si/Si02 interface results in the 
electron wavefunction penetration through the interface 
and into the oxide. The wavefunction penetration is ac- 
counted for by including a few mesh points in the ox- 
ide while solving the Schrodinger equation. ARPACK 
package^ was used to solve the 2D Schrodinger equa- 
tion and the successive over-relaxation (SOR) method 
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was used to solve the 2D Poisson equation. The conver- 
gence of the coupled Schrodinger-Poisson solver is found 
to be faster when the Poisson equation is solved by using 
the SOR method than when the incomplete lower-upper 
(ILU) decomposition method is used. 

The Monte Carlo transport kernel is used to simulate 
the electron transport along the axis of the wire under 
the influence of the confining potential in the transverse 
directions and a very small lateral electric field along the 
channel. The long wire approximation implies that the 
transport is diffusive (the length exceeds the carrier mean 
free path), and therefore justifies the use of the Monte 
Carlo metho d 42 ! 43 to simulate electron transport. Elec- 
trons are initialized such that their average kinetic energy 
is (1/2)KbT (thermal energy for ID) and are distributed 
among different subbands in accordance with the equilib- 
rium distribution of the states obtained from the Poisson- 
Schrodinger solver. Since the electrons are confined in 
two transverse directions, they are only scattered in ei- 
ther the forward or the backward direction; consequently, 
just the carrier momentum along the length of the wire 
needs to be updated after each scattering event. Mobility 
is calculated from the ensemble average of the electron 
velocities.— 



B. Scattering due to Bulk Acoustic Phonons, 
Intervalley Phonons, and Surface Roughness 

Phonon scattering and the SRS are considered in this 
work. The SRS was modeled using Ando's model^ in- 
tervalley scattering was calculated using bulk phonon ap- 
proximation, and the intravalley acoustic phonons were 
treated in both the bulk-mode and confined-mode ap- 
proximations. Since the wire is very lightly doped, the 
effect of impurity scattering was not included. Non- 
parabolic band model for silicon, with the nonparabol- 
icity factor a = 0.5eV _1 , was used in the calculation of 
scattering rates. A detailed derivation of the ID scatter- 
ing rates is given in Appendices S (phonon scattering) 
andlBl (SRS). Here, for brevity, only the final expressions 
for the scattering rates are given. 

For an electron with an initial lateral wavevector k x 
and parabolic kinetic energy £k . = S 2 ^/(2m*) in sub- 
band n [with subband energy £ n and electron wavefunc- 
tion ij) n {y, z)\, scattered to subband m [with subband en- 
ergy £ m and electron wavefunction ip m (y, z)], the final 
kinetic energy £f is given by 



£ r — £„ 



£ m H „ 1- nuj, 

la 



(1) 



where Tim = for elastic (bulk intravalley acoustic 
phonon and surface roughness) scattering, hco = ±Tkoq 
for the absorption/emission of an approximately disper- 
sionless intervalley phonon of energy TiuiQ, while in the 
case of confined acoustic phonons (below) the full phonon 
subband dispersion is incorporated. 



The intravalley acoustic phonon scattering rate due 
to bulk acoustic phonons is given by 



Z 2 ac k B TV2r 

h 2 pv 2 



(1 + 2a£ f ) 
y/£f(l + a£ f ) 



(2) 

where s ac is the acoustic deformation potential, p is the 
crystal density, v is the sound velocity, and 9 is the 
Heaviside step-function. V nm represents the overlap inte- 
gral associated with the electron-phonon interaction (the 
so-called electron-phonon wavefunction integral^) , and is 
given by 



V — 



\ip n (y, z) | 2 1 ip m (y, z)\ 2 dy dz. 



(3) 



The intervalley phonon scattering (mediated by short 
wavelength acoustic and optical phonons) rate is given by 



\f2T1pujQ 



N + - qp - 
2^2 



Vnr. 



(1 + 2a£ f ) 
y/£ f (l + a£ f ) 



(4) 



G(£ f ), 



where Sj„ is the intervalley deformation potential, and 
T> nm is defined in ([3]). The approximation of dispersion- 
less bulk phonons of energy hivo was adopted to describe 
an average phonon with wavevector near the edge of the 
Brillouin zone and Nq = [exp(hujo / UbT) — 1] 1 is their 
average number at temperature T. 

Assuming exponentially correlated surface 
roughness 4 ^ and incorporating the electron wave- 
function deformation due to the interface roughness 
using Ando's model^ the unscreened SRS rate is given 
by 



A 2 A 



2+( g ±) 2 A 2 
(l + 2ag/) 
\J£ f (l + a£ f ) 



(5) 



Q(£ f ), 



where A and A are the r.m.s. height and the correlation 
length of the fluctuations at the Si-Si02 interface, respec- 
tively, q^r = k x ±k' x is the difference between the initial 
(k x ) and the final (k' x ) electron wavevectors and the top 
(bottom) sign is for backward (forward) scattering. The 
SRS overlap integral in Eq. (jSJ) due to the top interface 
for a silicon body thickness of t y is given by 



nm. — 



dy dz 



&ty TYliy 



^ m (y,z) 



d 2 ip n {y,z) 



+ *Pn(y,z)e y (y,z) 

f £ - £ 

+ Vn{y,z) 



ty 



1 



dy 2 

ipm(y,z) 

y \ d<ip m (y,z) 



(6) 



dy 



The SRS overlap integral was derived assuming the 
interfaces to be uncorrelated. For the bottom interface, 
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the integration should be performed from the bottom in- 
terface to the top interface and the integral for the side 
interfaces can be obtained by interchanging y and z in 
Eq. ©. The first term in Eq. © is the confinement- 
induced part of the SRS and it increases with decreasing 
wire cross section. This term does not depend on the 
position of the electrons in the channel and hence results 
in high SRS even at low transverse fields from the gate. 
The second and third terms in Eq. ((6|) depend on the 
average distance of electrons from the interface, so they 
contribute to the SRS only at high transverse fields from 
the gate. (In Sec. IIII1 for example, we will see that the 
major contribution to the SRS in wires thinner than 5 x 
5 nm 2 comes from the confinement-induced term in Eq. 
(J6j> , and it increases rapidly with decreasing wire cross 
section.) 

Scattering rates given by Eqs. ©-((S]) are calculated 
using the wavefunctions and potential obtained from the 
self-consistent Poisson-Schrodinger solver. For this work, 
the parameters used for calculating the intervalley scat- 
tering were taken from Ref. 13a. the acoustic deformation 
potential was taken from Ref. [3fJ, and A = 0.3 nm and A 
= 2.5 nm were used to characterize the SRS due to each 
of the four interfaces. The SRS parameters were obtained 
by fitting the mobility of an 8 x 8 nm 2 SiNW in the high 
transverse field region (where the SRS dominates) with 
the corresponding mobility observed in ultra-thin SOI of 
similar thickness^ 



C. Acoustic Phonon Confinement 



In ultra-thin and ultra-narrow structures, the acous- 
tic phonon spectrum is modified due to a mismatch of the 
sound velocities and mass densities between the active 
layer and the surrounding material ) 20 ! 45 in our case - sil- 
icon and SiC>2. This modification in the acoustic phonon 
spectrum becomes more pronounced as the dimensions of 
the active layer become smaller than the phonon mean 
free path, which is around 300 nm in silicon4£ Pokatilov 
et al£Q have shown that the modification in the acoustic 
phonon dispersion in nanowires can be characterized by 
acoustic impedance C = pV s , where p and V s are the mass 
density and sound velocity in the material, respectively. 
By considering materials with different Q, Pokatilov et 
al£Q have shown that the acoustic phonon group veloc- 
ity in the active layer is reduced when an acoustically 
soft (smaller £) material surrounds an active layer made 
of acoustically hard (higher £) material. Since Si is acous- 
tically harder than SiC>2, the acoustic phonon group ve- 
locity in SiNWs with SiC>2 barriers decreases and results 
in an increased acoustic phonon scattering rate [see Eq. 
©]■ 




Longitudinal phonon wavevector [x 10 6 cm 1 ] 

FIG. 2. (Color online) Confined acoustic phonon dispersion 
(dilatational mode) calculated using the xyz algorithm^ 2 - for 
an 8 x 8 nm 2 SiNW. Only the lowest 10 phononic subbands 
are shown. Dispersion in the first one third of the first Bril- 
louin zone is shown for clarity. 

1. Confined Acoustic Phonon Dispersion 

The first step in accounting for the acoustic phonon 
confinement in mobility calculation is to calculate the 
modified acoustic phonon dispersion. Using the adia- 
batic bond charge model*£ (microscopic calculation, ac- 
curate but computationally involved), Hepplestone and 
Srivastava have shown the validity of the elastic contin- 
uum model (macroscopic calculation, less accurate but 
easier to implement) for wire dimensions greater than 
2.5 nm4£ Hence, in this work we have used the elastic 
continuum model to calculate the modified phonon spec- 
trum. Most of the previous studies of acoustic phonon 
confinement in nanowires have used approximate hybrid 
modes proposed by Morse 4 ^ (valid for wires with thick- 
ness much smaller than the width) to calculate the dis- 
persion spectrum. Nishiguchi et al~2- calculated the dis- 
persion spectrum using the xyz algorithm^! and found 
that the Morse formalism is valid only for the lowest 
phonon subband. Since one acoustic phonon subband 
is certainly not enough to accurately describe scattering 
with electrons, in this work, we have used Nishiguchi et 
aZ.'s approach to calculate the acoustic phonon disper- 
sion, although it is computationally intensive. The basis 
functions used to expand the phonon mode displacements 
in Nishiguchi et al.'s approach are powers of Cartesian co- 
ordinates, and the number of basis functions required to 
fully describe the modes depends on the number of modes 
required. For an 8 x 8 nm 2 SiNW, we have found that 
the lowest 35 phononic subbands are enough to calculate 
the scattering rate. Also, we found that about 176 basis 
functions are sufficient to fully describe the displacement 
of those 35 phononic modes. The number of phononic 
bands required decreases with a decrease in the wire cross 
section. 

Two types of boundary conditions are often used 
to calculate the acoustic phonon spectrum in nanostruc- 
tures: a) The free-standing boundary condition (FSBC) 
assumes that all the surfaces are free, so normal com- 
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ponents of the stress tensor vanish at the surfaces, and 
b) The clamped-surface boundary condition (CSBC) as- 
sumes that the surfaces are rigidly fixed, so the displace- 
ment of phonon modes is zero at the surfaces. Generally, 
the CSBC (FSBC) results in higher (lower) phonon group 
velocity than the bulk case^ For the wires considered 
in this work, neither of these boundary conditions holds 
exactly, since these wires are actually embedded in the 
Si02- Ideally, one needs to solve the elastic continuum 
equation, taking into account the continuity of displace- 
ment and stress at all Si-Si02 interfaces, and then apply 
the boundary conditions at the outer surfaces. But, this 
is almost numerically impossible for the structure con- 
sidered because it would be equivalent to solving the 2D 
Schrodinger equation in a device with the cross section 
of about 800 x 400 nm 2 , three (five if the metal-Si in- 
terface is included) interfaces along the depth, and two 
interfaces along the width. Donetti et aZ.)2S in their work 
on a Si02/Si/SiC>2 sandwich structure, considered con- 
tinuity of the displacement and stress at the interfaces 
to calculate the phonon dispersion. They found it to 
be close to the results from the FSBC. Therfore, in this 
work, we have used the FSBC to calculate the acoustic 
phonon spectrum of SiNWs. Fig. [2] shows the calculated 
acoustic phonon dispersion of the lowest 10 dilatational 
modes for an 8 x 8 nm 2 SiNW. Apart from these di- 
latational modes, depending on the rotational symmetry 
of the confined acoustic phonon displacement, there are 
two sets of flexural modes and one set of torsional modes 
in SiNWs. A detailed description of the symmetry of all 
these phonon modes can be found in Ref. [321 . 

2. Scattering due to Confined Acoustic Phonons 

The modification of the acoustic phonon dispersion 
due to confinement, shown in Fig. [2J implies that the lin- 
ear dispersion and elastic scattering approximation can 
no longer be used in calculating the scattering rate. The 
modified scattering rate which takes into account con- 
fined acoustic phonon modes is given by 

J i=l,2 v 7 



where q x is the lateral wavevector of the acoustic phonon, 
g(q x ) = E — E' ^-Tiuj j {q x ) , q Xl and q X2 are the two possible 
roots of g(q x ) = 0, and g'(q Xl ) and g'(q X2 ) are the deriva- 
tives of g{q x ) with respect to q x evaluated at q Xl and 
q X2 , respectively. Index J stands for the different acous- 
tic phonon modes and Nj qa . is the number of acoustic 
phonons of energy hu>j(q x ). Overlap integral C nm (J,q x ) 
and the total energy of the electron before (E) and after 
(E') scattering are defined in Appendix I A 31 

For intrasubband transitions, only dilatational 
modes are important, because for all other modes the 
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FIG. 3. (Color online) Electron - acoustic phonon scatter- 
ing rate for the lowest electron subband of an 8 x 8 nm 2 SiNW 
at the channel sheet density of N s = 8.1 x 10 11 cm -2 , calcu- 
lated assuming the bulk (dash-dot green line) and confined 
phonons (solid red line). The electron- bulk acoustic phonon 
intersubband spikes are at around 20 meV, 52 meV, 85 meV, 
and 95 meV, and they correspond to the electron scattering 
from the lowest subband to the 2 nd , 3 rd , 4 th , and the 5 th sub- 
bands, respectively. To the left /right of each intersubband 
scattering spike that corresponds to the bulk phonon approx- 
imation (dash-dot green) are two groups of small spikes (solid 
red) that correspond to absorption ("Ab") /emission ("Em") 
of confined phonons from different phonon subbands. 
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FIG. 4. (Color online) Group velocity of dilatational modes 
for an 8x8 nm 2 SiNW. On average, acoustic phonon group 
velocity is reduced to less than 50% of the bulk value (9.13 
xlO 5 cm/s). 

overlap integral C nm (J,q x ) in Eq.([7]) vanishes due to 
symmetry. In intersubband transitions, all the four sets 
of acoustic phonon modes are included in the calcula- 
tion of the electron - confined acoustic phonon scatter- 
ing rates, but the dominant contribution to the scattering 
rate comes from the dilatational modes. 

Fig. [3] shows the intrasubband electron-acoustic 
phonon scattering rate for the lowest electron subband, 
calculated using both the bulk-mode and confined-mode 
approximations. When calculating the electron-bulk 
acoustic phonon scattering rates, acoustic phonon disper- 
sion is assumed to be linear, w q — V s q, where V s is the 
sound velocity, as before. The resulting scattering rate, 
in the elastic and equipartition approximations (see Ap- 
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pendix Al), is proportional to the final electron density 
of states, and has the characteristic f D density-of-states 
peaks (dot-dashed green line) whenever the electron en- 
ergy becomes sufficient to scatter into the next subband. 
In the case of confined acoustic phonons, as seen in Fig. 
[21 the elastic approximation for electron-phonon scatter- 
ing no longer holds, and neither does the linear disper- 
sion at small wavevectors (except for the lowest phonon 
subband). Still, one can speak of a group velocity asso- 
ciated with a collection of phononic subbands. The aver- 
age group velocity accounting for the non-uniform energy 
gap between different phonon modes is shown in Fig. 2] 
The average group velocity is close to the bulk value for 
very small phonon energies, but asymptotically reaches 
a constant value (less than 50% of its bulk value) at high 
phonon energies. Since the scattering rate due to con- 
fined phonon subbands is inversely proportional to their 
group velocity, on average, the confined acoustic phonon 
scattering rate is about two times the acoustic scattering 
rate calculated using bulk phonons (Fig. [3]). Moreover, 
each of the bulk-phonon-scattering intersubband peaks in 
Fig. [3J (obtained in the elastic approximation, so the ab- 
sorption and emission rate peaks coincide) splits into two 
groups of peaks when confined phonons are considered: 
confined phonons can generally not be treated as elastic, 
hence there is a group of small peaks due to confincd- 
phonon absorption below each bulk-phonon- scattering 
peak associated with a given electron subband, and a 
group of peaks due to confined-phonon emission above 
the bulk-phonon-scattering peak. 



III. ELECTRON MOBILITY IN SILICON 
NANOWIRES - SIMULATION RESULTS 

In our previous work^ we examined the variation 
of the transverse field-dependent mobility with decreas- 
ing channel width in 8 nm thick rectangular SiNW of 
different widths and observed two that (i) the mobility 
at low-to-moderate transverse fields, limited by phonon 
scattering, decreases with decreasing wire width, and (ii) 
the mobility at high transverse fields, which is dominated 
by the SRS, increases with decreasing wire width. The 
former is due to the increase in the electron-phonon wave- 
function overlap [Eq. ([3])] with decreasing wire width, 
and the latter is due to the onset of volume inversion. In 
this section, we first emphasize the importance of acous- 
tic phonon confinement in SiNWs, and then vary the 
cross section of the wire to investigate the effect of in- 
creasing spatial confinement on electron mobility. 

The electron mobility in an 8 x 8 nm 2 wire, with and 
without phonon confinement, is shown in Fig. [5] In the 
low-transverse-field region, the mobility calculated with 
confined acoustic phonons is about 10 % lower than that 
obtained with bulk phonons. This clearly indicates that 
confined acoustic phonons need to be properly included 
in the study of electrical transport in SiNWs. The mo- 
bility values for 8x8 nm 2 are very close to the exper- 



imentally observed mobility in ultra-thin SOI of similar 
thickness^ In the remainder of the article, we will al- 
ways assume confined acoustic phonons. 
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FIG. 5. (Color online) Variation of the field-dependent mo- 
bility for an 8 x 8 nm 2 SiNW assuming bulk acoustic phonons 
(solid line) and confined acoustic phonons (dashed line). 

To determine the cross sectional dependence of the 
electron mobility in SiNWs and to understand the con- 
finement effects on the spatial and fc-space distribution 
of electrons, the cross section of the wire was varied from 
8x8 nm 2 to 3 x 3 nm 2 . The variations of the elec- 
tron mobility with decreasing wire cross section at a low 
(1.4 x!0~ 2 MV/cm), moderate (2.4 xIO" 1 MV/cm) and 
high (1.04 MV/cm) transverse field are plotted in Fig. [6] 
In the following, we will see that a complex interplay of 
several competing physical mechanisms is responsible for 
the electron mobility behavior observed in Fig. [5J 
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FIG. 6. (Color online) Variation of the electron mobility 
with SiNW cross section at three different transverse fields. 



A. Subband Modulation 

One of the most important factors that determines 
the energy and occupation probability of a subband in 
each of the Aq valleys (equivalent in bulk silicon) is the 
effective mass in the direction of confinement. For the 
SiNWs considered, the confinement is along the y and z 
directions and the electrons are allowed to move freely 
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in the x direction; consequently, the conductivity effec- 
tive mass for the valley pairs with minima on the x, y, 
and z axes are mi,mt, and m tl respectively, while their 
subband energies are roughly proportional to l/y/ mf , 
1/y/mtmi, and l/y/mtmi, respectively. Since rat < mi, 
the subbands in the valley pair along x are higher in en- 
ergy than those in the valley pairs along y and z. So the 
subbands split into those originating from the twofold 
degenerate A2 (the valley pair along x) and those origi- 
nating from the fourfold degenerate A4 valleys (the valley 
pairs along y and z). 

Upon increasing spatial confinement by decreasing 
the wire cross section, the subbands in different valleys 
are pushed higher up in energy, and consequently only a 
few of the lowest subbands in each of the valley pairs get 
populated with electrons. Fig. [7] shows the depopulation 
of the higher A2 valley subbands with increasing spatial 
confinement: since the lowest subbands in the A4 valleys 
are lower in energy than those in the A2 valleys, under 
extreme confinement A2 subbands get completely depop- 
ulated, and only the lowest A4 subbands are populated. 
Splitting of the valley degeneracy and modification of the 
subband energies in different valley pairs due to spatial 
confinement, followed by depopulation of the higher sub- 
bands, are together termed subband modulation. 35 Sub- 
band modulation enhances the electron mobility becauses 
it suppresses intersubband and intervalley scattering, as 
shown for ultrathin-body SOI MOSFETs^l In our pre- 
vious work on SiNWs^ we also observed a small en- 
hancement in mobility for wires of cross section around 
4x4 nm 2 , but in that study we did not consider the 
confinement-induced term in SRS. However, this term in 
the SRS increases rapidly for wire cross section below 5 x 
5 nm 2 and suppresses the beneficial effect that subband 
modulation has on the electron mobility in ultra-small 
SiNWs. Indeed, the confinement-induced SRS term has 
been shown to be dominant in determining the SRS in 
small cylindrical SiNWs^ 



B. Volume Inversion 

As the cross section of the SiNW decreases, the chan- 
nel electrons are distributed throughout the silicon vol- 
ume as opposed to just within a thin channel at the 
Si/SiC>2 interface right below the gate, as with conven- 
tional MOSFETs. The transition from surface inversion 
to volume inversion occurs gradually and the cross sec- 
tion at which the entire silicon is inverted depends on 
the electron sheet density— Fig. [8] shows the variation 
of the electron density across the wire with varying wire 
dimensions (gate is on top). When the cross section is 
decreased from 8x8 nm 2 (bottom right panel) to about 
6x6 nm 2 (bottom left panel), the onset of volume inver- 
sion results in an increase in the average distance of the 
electrons from the top interface, where the electric field 
is highest, so SRS is reduced [the second and third terms 
in the SRS overlap integral (|6|) drop with a decrease of 
the avarage electronic position from the interface, and 
are therefore sensitive to volume inversion]. But, once 
the wire cross section reaches about 6x6 nm 2 , the sili- 
con volume is fully inverted, so further reduction of the 
cross section simply results in a decrease of the average 
distance of the electrons from the interfaces (all four of 
them), thereby resulting in more surface- roughness scat- 
tering. Consequently, for wires with the cross section 
smaller than 6x6 nm 2 , volume inversion does not offer 
an advantage to electronic transport. 
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FIG. 7. (Color online) Variation of the electron population 
at a low transverse field (1.4 xlO -2 MV/cm). The solid line 
shows the total population of electrons in the lowest subband 
in the A4 valleys and the dashed line shows population of the 
A2 valley pair with varying spatial confinement. 



FIG. 8. (Color online) Electron density across the nanowirc 
at a high transverse field (1.04 MV/cm). When the cross 
section is reduced from 8x8 nm 2 (bottom right panel) to 3 
x 3 nm 2 (top left), the onset of volume inversion is evident. 
The color bar on the right is in 5 xlO 18 cm -3 . 



C. Mobility Variation with the SiNW Cross 
Section — the Big Picture 

Fig. [6] shows the variation of the electron mobility 
when the SiNW cross section is varied from 8x8 nm 2 
to 3 x 3 nm 2 , for a high (red), moderate (blue), and low 
(green) transverse electric field from the gate. 

Low and moderate transverse fields: with decreas- 
ing wire cross section, the intrasubband phonon scatter- 
ing increases due to the increase in the electron-phonon 
overlap integral (Fig. [5]); intersubband scattering and 
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intervalley phonon scattering decrease due to subband 
modulation; SRS increases due to the increase in the 
first term in the SRS overlap integral with increasing 
confinement. Overall, the mobility decreases with de- 
creasing wire cross section, very weakly for larger wires 
and much more rapidly for wires roughly smaller than 
about 5x5 nm 2 . As the wire cross section increases 
above above 7x7 nm 2 , we observe a very weak mobil- 
ity variation that results from the competition between 
an increase in intersubband scattering (number of occu- 
pied subbands increases) and a decrease in intrasubband 
scattering (electron-phonon overlap integral decreases). 
A similar weak dependence of the electron mobility with 
increasing cross section has been reported by Jin et al<^ 
for cylindrical SiNWs with diameters greater than 6 nm. 

High transverse fields: As the wire cross section is 
reduced from 8x8 nm 2 to 5 x 5 nm 2 , the first term 
in the SRS overlap integral ^ increases, whereas the 
second and the third terms decrease due to the onset of 
volume inversion (Fig. |HJ). Consequently, the mobility 
shows a very small change for these cross sections. But, 
when the wire cross section is smaller than 5x5 nm 2 , 
the benefits of volume inversion are lost and all the terms 
in the SRS overlap integral increase with decreasing wire 
cross section. 

Transverse-field independence of the electron mo- 
bility for very thin wires: We also notice that the 
transverse-field dependence of the electron mobility 
weakens with decreasing wire cross section and becomes 
virtually unimportant for SiNWs thinner than about 5 x 
5 nm 2 . Irrespective of the effective field from the gate, 
the mobility decreases monotonically with increasing spa- 
tial confinement, with the limiting mechanisms being the 
steady increase in the field-independent, confinement- 
induced part of the SRS © and the increase in intra- 
subband phonon scattering. Of these two mechanisms 
that limit the electron mobility in ultra-thin wires, the 
SRS scattering dominates. 
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FIG. 9. (Color online) Variation of the electron-phonon in- 
trasubband scattering overlap with SiNW cross section at a 
low transverse field (1.4 xlO -2 MV/cm). This results in a 
fourfold increase in the intrasubband phonon scattering rates, 
as the thickness and width of square wires are varied from 8 
nm to 3 nm. 



IV. CONCLUSION 



In ultrathin SiNWs, both electrons and acoustic 
phonons experience 2D confinement. In wires surrounded 
by Si02, an acoustically softer material, the acoustic 
phonon group velocity is lowered to almost half of its bulk 
silicon value, and leads to enhanced electron-acoustic 
phonon scattering rates. The electron mobility calcu- 
lated while accounting for the modification to the acous- 
tic phonon spectrum due to confinement is about 10% 
lower than the mobility calculated with bulk acoustic 
phonons. This result clearly emphasizes the need to ac- 
count for the acoustic phonon confinement when calcu- 
lating the electrical properties of SiNWs. 

We systematically account for confined acoustic 
phonons in the calculation and find that the mobil- 
ity decreases with decreasing wire cross section, very 
weakly for thicker wires and much more rapidly for wires 
roughly thinner than about 5x5 nm 2 . Also, we find that 
the transverse-field dependence of the electron mobility 
weakens with decreasing wire cross section and becomes 
virtually unimportant for SiNWs thinner than about 5 x 
5 nm 2 . 

For thicker wires at low and moderate transverse 
fields, the slow decrease of mobility with decreasing wire 
cross section is governed primarily by the interplay be- 
tween the beneficial subband modulation (less intersub- 
band and intervalley scattering) and the detrimental in- 
crease in intrasubband scattering. At higher fields, how- 
ever, the weak mobility variation with decreasing cross 
section stems from the competing influences of the benefi- 
cial volume inversion and the detrimental enhancement of 
the confined-induced surface-roughness scattering term. 

For very thin wires (below the 5x 5 nm 2 cross sec- 
tion) , the mobility decreases monotonically with increas- 
ing spatial confinement and becomes virtually indepen- 
dent of the transverse electric field. This occurs primarily 
due to the increase in the field- independent, confinement- 
induced part of the SRS © , and secondly due to the in- 
crease in intrasubband phonon scattering. In contrast to 
bulk MOSFETs, in which the surface roughness scatter- 
ing plays an important role only for high fields from the 
gate, electrons in very thin SiNWs are strongly influenced 
by the roughness regardless of the transverse field. This 
finding is important both for FETs with multiple gates, 
such as the FinFET,— as well as for ungated ultrathin 
wires used for thermoelectric applications^ or intercon- 
nects. 

This work has been supported by the National Sci- 
ence Foundation through the University of Wisconsin 
MRSEC. 
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APPENDIX A: PHONON SCATTERING 
1. Bulk Acoustic Phonon Scattering Rate 

The electron wavefunction in subband n, taking into 
account confinement along the y and z directions and free 
motion along x, is 



(Al) 



The displacement field u due to longitudinal phonons 
in the second-quantization form can be written as 



u(r) = £ 



2pQuj q 



(a a e lqr + ale- iqr ) 



(A2) 



where, a q and a q are the phonon annihilation and cre- 
ation operators respectively, Q is the volume, p is the 
density, and e q is the polarization vector. 

The perturbing potential which goes into the matrix 
element calculation is given by 



H ac = ^acV • u, 



(A3) 



where S ac is the acoustic deformation potential and u is 
the phonon displacement given by Eq. (IA2p . 

From Eqs. (|A2|) and (|A3|) . we can write the perturb- 
ing potential as 



H QC = E ac £ y ^^e q • zq (a q e^ - a^) . (A4) 

The matrix element after integrating over phonon coor- 
dinates, for scattering from subband n with wavevector 
ka; to subband m with h x is given by 



M nm (k 



x- ,e 



My^y [qyV+q ' z H m {y,z)\ dydz ( A5 ) 



where 2V q is the number of phonons given by the Bose- 
Einstein distribution function 



vv q = 



(A6) 



Defining l nm {q y ,q z ) as 

I nm (q y ,qz) = JJ [ip n (y,zy {qvV+q * z) i>m(y,z)] dy dz 

(A7) 



and after integrating over x, Eq. (|A5[) yields 
M nm (k ;c , k^,) 



2 ft 



'2pQ,uj^ V " 2 ' 2 J (A8) 
x \l nm (q y ,q z )\ 2 S(k x - k x ±q x ). 

In the equipartition approximation, the phonon number 
becomes 

K B T 



(A9) 



so using the equipartition approximation and linear dis- 
persion relation for acoustic phonons defined by w q — v s q 
in Eq. (|A!j) . we get 



„ 2 hq 2 K B T 
2pflv s q hv s q 



2E 



(A10) 



x \l nm (q y ,q z )\ S(k x -k x + q x ), 



where v s is the sound velocity in the crystal. 

Acoustic phonon scattering rate using Fermi's 
Golden Rule with the elastic scattering approximation 
is given by 



IT m (kx) \M nm (k x ,k x ) '6(S-8' 



(All) 



9||>fc« 



where £ and £' are the initial and final energies of the 
scattered electron in the parabolic band approximation, 
respectively. 

Substituting Eq. (jAlOjl in Eq. (jAllj) and changing 
the sum to an integral, we get 



r^mAkx) 



2Tittz 2 „K B T ff ,:2dq y dq z 
\l nm {q y ,q z )\ 4?r2 



pQv 2 h 

J S (k x -k x + q x )s(£-£')^- 



To evaluate Eq. (|A12[) . let us rewrite it as 



2nZ a c 2 K B T 
hpv 2 s 



where 



J 2 = — / 8{k x - k x + q x )5{£ - £ ) dk x ', 
Zir 



Vnrn = l^nm {q y , q Z )\ 2 dq y dq z . 

The above equation can be written as 
1 



(A12) 
(A13) 

(A14) 
(A15) 



'At: 2 



dq y / dq 



xjj n (y,z)e iiqvV+q * z) iPm(y,z)] dydz 

My,z')e-i w ' +qJ )^ m (y, z ') 



dy dz . 
(A16) 
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Using the identity 

±-jdq y e^y-y^=S(y-y), (A17) 
we can write (|A16[) as 

P,„„ I I \ip n (y,z)\ 2 \^ m (y,z)\ 2 dydz. (A18) 



Adding a nonparabolicity factor a and converting 
the dk x integration to d£' integration, the integral in Eq. 
(IA14|) simplifies to 



1 + 2a£ f 



2vr V 2h 2 y/£ f (l + a£ f y 



(A19) 



where £f is the final kinetic energy of the electron af- 
ter scattering. It is defined in terms of initial parabolic 
kinetic energy £ kx as 



£f — £n £m 



VI + Aa£ kx - 1 
2a 



(A20) 



Substituting the above simplified integrals in Eq. 
(fAT3l) we get 



where the emission and absorption of an optical phonon 
results in (A^ q + 1) 5 and (N^)^ , respectively, and 
Inm (<ly , <lz) is the overlap integral defined in Eq. I|A7[) . 

Following the procedure outlined above to calculate 
acoustic phonon scattering rate and accounting for the 
inelastic nature of scattering due to optical phonons, we 
can write the optical phonon scattering rate as 



M = 



1 1 

2 T 2 



(ly><lz)\ 



• dq y dq z 
4tt 2 



8 (k x - k x ±q x ^j S(£ 



(A25) 

where S(£ — £ ±huj ) ensures the conservation of energy 
after absorption (top sign) and emission (bottom sign) of 
a phonon of energy Tiloq. Simplifying this using the ap- 
proach followed in the previous section (|A 1[) . intervalley 
phonon scattering rate can be written as 



_ ^ a c 2 K B TV^ 
l nm (k K j- ^_ U n 



1 + 2a£ f 



,®{£ f ). 



^£ f (l + a£ f ) 

where Q(£f) is the Heaviside step function that ensures 
a positive kinetic energy after scattering. 



2. Intervalley Scattering Rate 



-pit; \ 



— 2 



2y/2%pu>o 



1 1 

2 T 2 



Vr, 



(1 + 2af/) 
y/£ f (l + a£ f ) 



(A26) 



where the final kinetic energy £/ is similar to that in Eq. 
(|A20p . with ±huo to account for absorption (top sign) 
and emission (bottom sign) of a phonon. 



Intervalley scattering can be mediated by long- 
wavevector acoustic phonons or non-polar optical 
phonons. Intervalley scattering is modeled using the non- 
polar optical phonon model. 

The perturbing potential which goes into the matrix 
element calculation is given by 

H iv (q) = E iv eq ■ u, (A22) 

where S„ is the intervalley deformation potential and u 
is the phonon displacement given by Eq. (|A2[) . 

Assuming intervalley phonons to be dispersionless, 
i.e., w q = luq, Eq. (|A22[) can be written as 



n iv = E iv E^/^K^ + ^n- (A23) 

The matrix element for intervalley scattering is given by 
M -( k - k ^=^\/5^(^ + rt) 2 (A24) 

x Inm (Qyi Qz) 3(k x ~ k x zkq x ), 



3. Confined Acoustic Phonon Scattering Rate 

Using the xyz algorithm, the normalized displace- 
ment components for the J th acoustic phonon mode in 
terms of a complete set of basis functions $\ can be writ- 
ten as 

uj,i = ajxj,i\^\, (A27) 

where i = (x,y,z) represents one of the components of 
the displacement, aj is the normalization constant and 
XJ,i\ are the coefficients of the basis functions. 

Taking the center of cross section of the wire as the 
origin, the basis functions in terms of powers of Cartesian 
coordinates in the lateral directions are 

*A(w)=(^y (A28) 

where A = (r, s) , q x is the longitudinal wavevector of 
the acoustic phonon mode along the axis of the wire, 
and W and H are the width and thickness of the wire, 
respectively. 
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Following the normalization procedure indicated in 
Ref. H2, we get the normalization constant in Eq. (|A27[) 
to be 



acj 



(A29) 



where L x is the length of the wire, uij is the frequency 
of the J th phonon mode, and \j is the eigenvector cor- 
responding to ujj. E is the matrix as defined in Ref. [32T 
The acoustic phonon field, which is used to determine 
the perturbing potential, is given by 



u = Et aj 9-~ +"';.. le, 



,t 



(A30) 



Considering Eq. (|A30|) instead of Eq. (|A2[) to repre- 
sent the phonon displacement, the matrix element given 
by Eq. (|A5[) can be rewritten as 

/ 1 1\5 

M nm (k x ,k x ) =^ ac a.j [Nj^ + \ C nm (J,q x ) 



1 



dx, 



(A31) 

where Nj qx — [exp(hu)j(q x )/kBT) — 1] is the number 
of acoustic phonons of energy Twoj(q x ) and C nm (J, q x ) is 
the electron-phonon overlap integral given by 



^nm^J-,qx) — 



2rx.J,i\ (lz\ r 1 ( 



W 



H 



w 

s-l 



w 



H 



>ip m {y,z) 



dy dz. 



The square of the matrix element is then given by 
2 



(A32) 



M nm (fe x , k X/ 



H* [ N. 



*Jq x 



1 1 

2 T 2 



x \C nm (J, q x )\ 5(k x ~ k x ±q x ). 

(A33) 

Considering the confined acoustic phonon scattering 
to be inelastic, the scattering rate can now be written as 



—2 

' — a 

2WH 



^E I ^[n^+It 1 

J 



2 2 



\<*J 



x \C nm (J, q x )\ S(E - E ±?ujj(q x )), 

(A35) 

where the total energy of the electron before (E) and 
after (E') scattering are defined as 



E = £„ 



E — £„ 



1 + 4a 



2m 



2a 



1 + 4a 



ft 2 (fc x ±g x ) 2 



(A36) 



2a 



The argument of the delta function in Eq. (|A35|) can 
have two roots. Using the identity for delta functions 
with multiple roots, the final expression for the scattering 
rate is 



—2 

2WH 



2 ' 2 



x \C nm (J,q x )Y 



8{q x - q x \) . S{q x - q x2 ) 



\g'{Qxi)\ \g'(Qx2)\ 

(A37) 

where g(q x ) = (E - E ^huj(qx)), q xl and qx2 ars the 
two possible roots of g{q x ) = 0, and g (q x i) and g (q X 2) 
are the derivatives of g(q x ) with respect q x evaluated at 
q x i and q X 2, respectively. 



APPENDIX B: SURFACE ROUGHNESS 
SCATTERING 

In a very simple model used to describe the surface 
roughness scattering, the perturbing potential for the in- 
terface normal to the y direction is given by 



H sr = ee y (y,z)A(x), 



(Bl) 



where A(x) is a random function which describes the de- 
viation of the actual interface from the ideal flat interface 
and e y (y, z) is the field normal to the interface. The scat- 
tering matrix calculated using this perturbing potential 
is 



^nm (^x) 



1 1 



^E(^ +:>>:>) K'l 

J,1x 



x \C nm (J,q x )\ S(k x - k x ±q x )8(E - E ±hujj(q x )), 

(A34) 

where the upper and lower signs denote absorption and 
emission of an acoustic phonon of energy ?iu>j(q x ), re- 
spectively. Integrating over q x and including the non- 
parabolicity factor, the scattering rate can be written as 



M nm (k x ,k x ) =e J J [ip n (y,z)e y (y,z)ip m (y,z)]dydz 
1 



x— / A(x)e 4(fe *- fc => dx. 



(B2) 

Defining the term in the double integral of the above 
equation as T nm , the square of the matrix element can 
be written as 
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2 _ 2 2 _J_ / i 

~ e •'nm £2 



x / dxA^A^V^'^ - ^- 



M ram (k a; , k Xj 



(B3) 

The average value of the matrix element given by Eq. 
(|B3|) over many samples is actually used to calculate the 
SRS. The expectation value of the square of the matrix 
element is given by 



M nm (k x , kjj.) 



J = e 2 T 2 — I ,/, 

c J nm jji ' 



x / dx(A{x)A(x'))e l( - x ~ x >•«-, 

(B4) 

where the correlation function R(x — x ) = (A(aj)A(a; )) 
depends only on the distance \x — x |, and q x — (k x — k x ). 
Redefining (x — x ) as x we get 



M lk,.k ([ ) ) = J_ / ,/., / ,/.,■""/?(.,•' )r" ■< 



(B5) 



Assuming exponentially correlated surface roughneS: 



,44 



defined by i?(x) — A 2 e * in Eq. (|B5|I . we get 



M nm (ka;, k^) 



— P 2 'T 2 A 2 T, 



2y/2X 



LI X {qf) 2 X 2 +2 



(B6) 

The scattering rate due to the interface imperfections 
can now be written as 



k'. 



L x q 2 X 2 + 2 



,(B7) 

Converting the sum over fc^ to integral over cZfc^, , we 



get 



1 



*2l j dK {k x ±k\^\- + 2 



(B8) 



<5(£ -£). 



Defining £, = ^ 1+4 °£>^ — and redefining (g*) 2 as 

2m/h 2 { v / £ i {l + a£ i ) ± ^(1 + a£/)) 2 , the final ex- 
pression for surface roughness scattering assuming non 
parabolic bands can be written as 



A 2 A 



2 + (g*) 2 A 2 
(1 + 2a£f) 



(B9) 



y/£f(l+0£f) 

Unlike the approach detailed above, Ando's models 
of interface roughness scattering accounts for deforma- 
tion of both the wavefunction and potential due to the 
imperfection at the Si-SiC>2 interfaces. The matrix ele- 
ment including the perturbation to wavefunction and the 
potential using Ando's model is given by 



M nm (k x ,k a .) = / / dydz I dx 



VI 



ip m (y + Ax, z) 



x [Ho + AV(y + Ax, z)] ip n (y + Ax, z) 
il> m (y,z)R ip n (y,z) 



(BIO) 

where the unperturbed system is represented by the 
Hamiltonian Ho and wavefunction ip n (y,z) and the per- 
turbed system's Hamiltonian and wavefunction are [Ho + 
AV(y + Ax,z)] and tp m (y+ Ax, z). Calculating the scat- 
tering rate from this matrix element, we find that the fi- 
nal expression for the SRS is same as before except that, 
we have additional wavefunction deformation terms in 
the SRS overlap integral T nrn . The SRS overlap integral 
is given in Eq. |6]). 
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